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Very accurate analytical expressions for the radiative corrections of unpolarized hyperons semilep- 
tonic decays of charged and neutral baryons have been obtained in the recent past. Some of these 
formulas contain logarithmic singularities at the edges of the Dalitz plot for the three and four-body 
decays. These singularities are analyzed and integrated analytically to obtain new divergentless 
formulas for the energy spectrum of the produced /3 particle. The new equations contain terms of 
, the order a times the momentum transfer, are applicable to any beta decay process and are suitable 

■ for a model-independent experimental analysis. 
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^ . The subject of hyperon decays is of unquestionable importance in high energy physics. The precise description 
of the semileptonic weak decays and in particular the neutron decay, is relevant in cosmology, astrophysics, solar 
physics, the solar neutrino problem, and in other areas of particle physics. A precise formula for the energy spectrum 
^-H ■ of the decay products (fcrmions) in the hyperon semileptonic decay (HSD) requires the knowledge of the radiative 
J> I corrections (RC) in the whole region of the Dalitz plot (DP). The total decay rate can be computed directly through 
CO . analytical formulas for the RC, if these formulas do not contain any divergences. The subject of this paper is to 
show how the logarithmic kinematical divergences, that appeared in the previously calculated model-independent 
bidimensional distributions in Refs. obtained in an analytical way for the RC to baryon p decays in the three 

body region (TBR) and the four body region (FBR) of the DP, are canceled after performing the integration over the 
, energy of the final hyperon that emerges in the process. The new results are valid in any charged or neutral hyperon 
beta decay. The important feature of the obtained analytical forms of the RC's is that they are written in a simple 
form as products of two factors. One of these factors is a model-independent function and the other one does depend 
f-^ ^ on the form factors whose determination (through the experimental data) is useful to obtain information about the 
I underlying interactions in the decay processes, basic symmetries and the internal structure of hadrons. The new 
O-i formulas are suitable for a direct evaluation of the RC's for any event at any point of the allowed physical region. 
^ The structure of this paper is as follows. In Sec. II we exhibit the kinematical region in which the three and the four 

body decays take place. All the amplitudes, the one without any RC, the bremsstrahlung and the virtual ones, for 
the charged and neutral HSD, with all the aq/irMi terms included, are displayed in this section. The bidimensional 
observable energy distributions in the FBR and TBR are described for the charged and neutral HSD processes in 
Sec. HI . We devote Sec. IV to display the formulas for the energy spectrum of the emitted charged lepton in a 
[ non-divergent form for the pure bremsstrahlung correction in the FBR, and the formulas for the energy spectrum 
in the TBR with both kinds of RC 's (the combination of the bremsstrahlung and the virtual contributions). In 
Sec. V we present our final formulas in each region for the charged and the neutral decay processes. Conclusions are 
given in Sec. VI, with the formulas that contain all the up-to-date improvements and new features for the energy 
spectrum, together with a comparison with other previously published numerical values assuming their values for the 
form factors. Several graphs for the partial radiative corrections (the Coulomb effect is included for the neutral HSD 
process at the TBR) and for the complete electron spectrum are presented in order to illustrate our results. Finally, 
for self-containment, we include several appendices with the definitions of the coefficients and model-independent 
functions that appear in the final analytical results, among other relevant relations. The procedure of integration 
with the explicit cancellations of the kinematical divergences is also included. 

II. KINEMATICS AND AMPLITUDES 

We start with the presentation of the main features and the notation to describe the decay processes we are 
interested in: 

A%p,)^B{p2)+e-{£) + DM, (1) 
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where the emission of a virtual or a real photon j{k) takes place. The emission of the photon is described as a 
radiative correction to the semileptonic decay of the hyperon. corresponds to the neutral (s = n) or charged 
(s = c) decaying baryon, B is the produced baryon, and denote the lepton and its antineutrino counterpart, 
respectively. The four-momenta and masses of the particles involved in the baryon semileptonic decay are denoted by 

Pi = {Ei,pi), P2^{E2,P2), (^(EJ), p^ = (E^,p^), andfc = (fco,fc), 



and by Mi, M2,m,m^, and , respectively. We assume throughout this paper that nii, = 0, and rn^ = for the 
real photons. 

The conservation of energy and momentum determines the physical region where the HSD takes place. For the 
three body decay as in the HSD (without any radiative correction) the physical region is bounded by the hyperon 
minimal and a maximal energy (^E™™, E^'^^^ for each value of E. When a real photon is considered as an additional 
product of the decay the physical four body region for this process contains the former region (TBR) and an additional 
portion (FBR) with energies below the £^™™(iJ) in such a way that M2 < E2 < i?™"' ior m < E < Ec where Ec is 
the maximum electron energy when the final _B-hyperon is at rest. The upper boundary of E2 , of the region (FBR) 
where the three body decay is forbidden, is evaluated with the condition fc = 0, this means, without the emission 
of a real photon. Physically this corresponds to the case in which the neutrino balances the total momentum of the 
residual hyperon and electron that are emitted in colinear directions. 

Consequently, the Dalitz Region for the process in Eq. in the rest frame of A^, is defined at the FBR by 



M2< E2< E. 



min 
2 I 



m< E < Er 



Er. 



(Ml 



2 (Ml - M2) 



(2) 



and at the TBR by 



£;™'" <E2< E'^''''., 



m < E < E„i, 



E„ 



Ml 



M| 



2Mi 



(3) 



where 



= -{Mi-E±pE) 



M' 



2 (Ml - E±(3Ey 



and 



^E. 



(4) 



The z-axis is chosen along the electron three-momentum and the x-axis oriented so that the final baryon three- 
momentum is in the first or fourth quadrants of the x — z plane. To be more explicit, 



Pi = (Mi,0), 



f = /?i?(0, 0, 1), 



P2 



m (V(i-y2), 0, y 



(5) 



In the TBR, the radiative correction is due to a virtual or a real photon emission, meanwhile in the FBR the emitted 
photon is a real photon exclusively, i.e., the radiative correction is due only to the brcmsstrahlung effect in the 
semileptonic decay. 

The uncorrected matrix element Mq (without any emission of a photon at the TBR) for the decay in Eq. (^ is 
given by the product of the matrix elements of the baryonic weak current and of the leptonic current: 



Mo = ^^UBW^j.UAUtOf^Vu 
V2 



(6) 



where = G^Vij and is the muon decay coupling constant, Vij is the corresponding Cabibbo-Kobayashi-Maskawa 
matrix element, and 



Ml). 
Ml 



Ml 

Ml 



gjil). 

Ml 



gjil). 

Ml 



-9m 



75, and Op =7^(1-1-75). (7) 



The q = pi — P2 denotes the four-momentum transfer. 

To obtain the precise decay rates, one has to consider all the involved amplitudes of the processes with and without 
RC. It has been shown in Ref. [Q that the order a amplitude can be obtained in a model independent fashion by 
using the Low-theorem, Refs. ||]. 

For completeness, we reproduce here the model independent brcmsstrahlung amplitudes, in terms of the Dirac form 
factors, given for the charged HSD in Eqs. (18)-(20) of Ref. 0, and for the neutral HSD in Eq. (23) of Ref. g. 
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The total hremsstrahlung transition amplitudes for the charged HSD (s = c) and the neutral HSD (s = n) decays 
in the entire Dalitz plot region are 



M% = Ml + + M|, with 



M2 = -^et^UBWxUAUi^^^^Ow,, 



(8) 



M3 = --j=ue,OxVyef,UB 
+e — ^ - (/^A 



2pi • k 2pi ■ k 

h - /2 , 53-52 

-75 



Ml 



Ml 



2p2 ■ fc 



/2 + .9275 \ 

i\fi j 



(9) 



and 



Ml" = eMo 



e • £ e • p2 
£ • A: P2 ■ k 



MJ,^ = M2^ 



M3 = —ueOxv„ef,UB < — ^ — 



\P2 ■ k 



9w 



/2 + .9275 
Ml 



'«iVVa— pcr^i/fc,. + K2<7^^k^— -PVa 

2pi ■ k 2p2 ■ k 

h + 9315 



<y\p + 9>.p 



Ml 



UA- 



(10) 



Ki and K2 are the anomalous magnetic moments of A* and B, respectively, is the photon polarization four-vector. 
The total wrtuaZ transition amplitudes (present only at the TBR) with terms up to order of aqj-KMi can be written 

as 



M^ = M^ + M,„ 



where Mq is the uncorrected matrix element in terms of effective form factors Ref. 
Explicitly, for the charged HSD case 



a 



zvr 



where 



arctanh (/?) — 1 



-— arctanh {(i) 
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2(3 



Ml - 2E 
1-/32 - 



- - In 



1 - 



2/3 



Mi/E- l + P 



In 


Ml 






m 



/3 



arctanh (/3) 



Ml - 2E 



(3E 



Ml - 2E 



:ln 



Mi/£; - l + /3_ 
arctanh (/?) 



Ml 



2^- 
11 



Ml 



with 



rx ^\ G E 

L(x) = / dt is the Spence function, and Mp. — ubW„uaui diOpV,^. 

Jo t V2 wMi 

The formulas for the neutral HSD are: 

= 7^[Mo$„ + Mp,$;], where Mp^ = -^—-p—-^UBWpUAUi^20pV^. 
2tt \/2m[l+p2Y 

The $„ and $^ are given in Ref. H, Eqs. (8)-(10), and can be written as 



(11) 



(12) 



(13) 



(14) 



(15) 



(16) 
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(17) 



where 



<^(£;,i?2,A) = 2 



'^niE,E2) 



where 



and 



-arctanh(/3jv) — 1 
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Pn Pn 
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-L 
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M2 



11 
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2A/| 



arctanh(/3Ar 
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1 



arctanh(/3Ar) 
M2 



In 



M| + (1 + pj,) a 

2{l+P2f 

— arctanh(/3jv) 



2M| 



Pn In 



M2 



/^AT = 1 ^ 



2A,f2x 1/2 



, a^2E{E2-P\p2\yo) = 2l.p2, 



- — - 



a/3. 



N 



2 /; I \2' 
+P2) 



+ i (1 ± /37v)a 

(^+P2)' 



(18) 



(19) 



(20) 



(21) 



(22) 



The contribution due to the Coulomb interaction in the neutral HSD is included in Eq. (|T^) . The infrared divergence 
which takes place in Eqs. (|l^) and (^8|) when A ^ in the virtual radiative correction is canceled by another one of 
the same kind which appears in the bremsstrahlung correction, see Ref. Q. A is the small mass which is assigned to 
the photon to cut off the contribution of very soft photons. 

In general, the ^s{E) and ^'s{E) are model independent functions (do not depend on the form factors) that contain 
terms of order aq/nMi. 



III. BIDIMENSIONAL DISTRIBUTIONS 

The evaluation of the differential decay rate that gives the DP with RC of process in Eq. (|l|) is performed by standard 
trace calculations, leaving as the relevant independent variables the energies E2 and E of the emitted baryon and 
the electron, respectively. The phase space integration is accomplished according to the kinematical limits given in 
Eqs. (|) and (|), for the TBR and PER respectively. 



A. Four body region 



In this region, the result up to order aq/irMi for the HSD is compactly given by products of coefficients iJf that 
depend on the form factors and model independent functions as 

dT's ( A' ^ Bei^ 7) dE = ^dQ ^ Hf Of , dVL = ^ —MidE2dE, (23) 

i=0 

where the coefficients Hf (see Appendix A) are given by 

HI ^ HI' = H[, H^ ^ Hi H^ =H[ + N[, for i = 0, 2, 17, H^.j = 0. (24) 

The 's do depend on the form factors through functions Qi, i — 1, ..,4. The explicit form of the Qi 's (Ref. [|o|) 
is also given in Appendix A. The model independent functions are: 
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QcT _ nriT 



for i = 0,2,...,17. 



(25) 



The ^J^, and 9j 's with i = 0, 17 are shown in the Appendix B. 

The specific form of Eq. (|2^ ) for the cases that correspond to the semileptonic decay of a charged charged HSD (s 
c) or a neutral neutral HSD (s = n) decaying hyperon, respectively, becomes: 
For the charged HSD process, such as S~(pi) n{p2) + e^(£) + Pe{Pv) + 7 (k), 



civ's ( ^ B°e-Pe 7) di; = ^ ij; 6*^ , 



16 



i=0 



where the model-independent function is 



9f 



ih - 2) In 



2/0 + 1 



2/0-1 



^1 = arctanh (/?) , ya 



2\p2\Ep 



, and E^ = Ml- E2- E. 



For the neutral HSD case, such as the A(pi) P^{P2) + e (^) + v^{pi,) + 7 (/c) 
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where 0"^ is separated, for further convenience, into two parts as follows: 



QuT _ nnTD , nnTND 



where 



'^nTD 



iE,E2 



-2 In 



yo + 1 



2/0-1 



(26) 



(27) 



(28) 



(29) 



The equations for the Q^'^^^ and 6*^7 are also included in the Appendix B. 

The result for both types of unpolarized decays given in Eq. (^3|) is of high precision, model independent and 
useful for processes where the momentum transfer is not small and therefore cannot be neglected. Terms of order 
aq^ /ttMI and higher are neglected. 



B. Three body region 



The result up to order aq/nMi for the hyperon semileptonic decay (HSD) in this region is given by 



d^TBR ( Bev ) dE = 



17 



Gl 1 



E2 



dn, dn^ — — ^Mi Af2 dEdz, —. (30) 



2 27r 



Ah 



The coefhcients 's are the same as the ones at the FBR (Appendix A), but the model independent functions are 
different, here 



Ai = B'{, Ai=AU, e[' = ^, + et, e[' = e[^ = e,, for z = o,2,...,i7. 



The ^0 and A^^. 's (given in Appendix A) also depend on the form factors through functions Qi, i — 1,..,4. The 
model independent functions O'l ^ , and the 6*1 's (i = 0, 2, 17 ) are exhibited in Appendix C. 

In order to illustrate the result in Eq. (|3^), let us consider two cases that correspond to the semileptonic decay of 
a charged and a neutral decaying hyperon: 
For the charged HSD 



dT\A- ^ B^e-V)dE 
For the neutral HSD, 



16 



A'^ + -\ Hl^eo + H[ ($e + ei) + B'l^', + E H[e, 



dfi. 



dr"(A° ^ B+e-V^) dE 
The 9[ " are analyzed in the Appendix C 
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dfi. 



(31) 



(32) 
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IV. ENERGY SPECTRUM 

To obtain the energy spectrum of the beta particle, the integration over E2 has to be performed in the corresponding 
kinematical regions, (i?26 < E2 < i?2t), according to Eq. (H) and Eq. (||) 



<'E2t 

The result for the case without any radiative corrections becomes 

-•2 ji.r rE\ 



K(A- ^ Bev) = |jf / Ai m. where |M„P = ^j^X^ (34) 

spin 



A. Pure Bremsstrahlung correction in the FBR 

The fir^(A* — > Bei^j) in Eq. ) depends on j/o through Of^ [see Eq. ( |27| ) and Eq. (p9|)] in the following way: 

2/0 + 1 



cx In 



(35) 



yo - 1 

For coUinear events, which take place when 

E2 Ef'', m<E<Ec (36) 

the HQ 1. This gives rise to a logarithmic divergence. This divergence is an obstacle to perform a direct numerical 
integration. However, this difficulty is controlled by performing an analytical integration. To proceed, it is convenient 
to take the following definitions: 

(yo + 1) ao + + s Ml {E^ -E) + M| E ^ Mi r-^ 

(yo - 1) ao + oqz - s M2E(3 E/3 

and to consider separately the terms where the difficulties arise. 



For the charged HSD the separation of terms in Eq. (26 ) is such that the ones that contain a logarithmic divergence 



(at the integration limits) are accumulated in and all the other non-divergent ones are accommodated in /^^^, 
then 

r%{E) = f [IE'"" (E) + (E)] , (38) 

where 

JTD ^ « f""' QCT ^ ^ ^c^jk j^TB^ (39) 

The explicit z dependence of H[ {E,E2) is given by 



H[=Y: 4 = ih - 2) s,, and = / In °° ^ + ] dz. (40) 



The Efc 's are shown in Appendix D. The second term in Eq. ( |38| ) is 
For the neutral HSD, we follow the same strategy and we obtain 



6 



where 



and 



(E) - - H[ef^ dE, = ^M, jz el M^Rr, el = -2s 



tTND 



(E) 



16 



A/2 



i=2 



Notice that, we obtain the same integral i?^^in both cases. After performing a very subtle analysis (illust 
Refs. Q and ^) we get the following non-divergent analytical result 



nTB _ 
* + 1 



rpr> 



r=0 



k — r 



dz 



where 



T,^ = In 



Afi [E.^ - E) 



EI3M2 



and Em is given in Eq. (||). The values of z and s (see Eq. ) at the boundaries are 



Zt 



St 



Mi-E{l- P) 

M2 

1 (Mi-E{l- (3) 



Mo 



Ml (1-/3) 



, Z, ^ Zb 



Ah-E{l + [3) 



Mo 



Mo 



2 V M2 Ml- E{1- (3) 

An analogous integral is discussed in Appendix E 



St = Sbi = 



M2 Ml -E{l + f3) 

Mi-E{l + P) M2 



Mo 



Ml -E{l + P) 



B. Virtual plus Bremsstrahlung corrections at the three body region 

Following a parallel development for the charged HSD we obtain 

r^(i?) = ^^Mi [/^ (E) + /^VD (E)] , 

where 



[E) = -M2 

TT 



.fe=0 



and 



^2 

I'nd{E)= J 



16 



+ - 1^(0^'' + B'^K + J^o^o + E ^^^^ 



i?fc(z) - J zHn 



ao 



+ hoz + v^^^ - 1 



dz. 
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The £j. 's are given in Appendix D, and the calculation of the integrated form of Rk{z) is illustrated in Appendix E. 
A similar procedure is performed for the neutral HSD, and the following results are obtained: 

r"(^) = ^^^^1 l^D (E) + i^D im , (51) 



(E) = -Ah 



(52) 



r r / 16 17 

Ind{e)= / j^, + '^[H[er' + A'i^'^'N + {H', + K)e, + Y,H[e, + Y,N[ 

^„i„ L \ j=2 z=2 



(53) 



The 9^^'^ and the other model independent functions are given in the Appendix C. 
Here we have considered the approximation 



arctanh(/3jv) I 



N 



di _ 1-1 ^ /^i _ !!f_!± \ - ll _ ^^200 /, _ m^h 
2 Mip \ 2) ^ 2 Ml/3 E^ 2 



Ml/? 



y 



(54) 



According to Eq. ( CIO ), the model dependent coefficients in Eq. {p2\j, defined as ej.", become 



J n _ J ^ 

£0 — — ^0 



1 



Ml/3 



£3 - £3 - 



Ml/3 



m 

1 ^/i 

2E^ 



£'i" = e'l + e'i = £1 ( /i - 2) - (5o£o + Af2aoei) 

= 4 + = £2 ( /i - 2) - (6o£i + M2ao£2) TF^ ( 1 - 



(55) 



The analytical result for the R^. 's in Eqs. (Esl) and (p^) is 



k+i ,fe+l^ 



2(fc + 1) 



-Tn 



2(fc + 1) 



fc-r+l _ A;-r+l\ 



(fc - r + 1) 



-dz 



(56) 



where 



In 



4 (£;„ - £;)' 



M| 



26 - Ec) In 



Ml (£;„ - E) 



M2EP 



(57) 



the 8 (i? — Ec) is the Heaviside function, and the zt, zf, St , and sf are given in Eq. (^q ) 



Zb = 



(58) 



V. FINAL RESULTS 
A. Four body region 

Gathering and refining previous results, we obtain the whole spectrum of events in the FBR. 
The FBR formula for the lepton energy spectrum becomes in general, 



8 



k=0 J^-^ \ i=2 ) 



dz 



(59) 



the (5" indicates that the last term appears only in the neutral HSD case (s = n). The ef. is given in Eq. ( |4C| ) for 
s = c, and in Eq. (^) for the s = n case. 



For /c = 0, 1, 2 the exphcit values for i?^^, in Eq. (^) are 



1 r 



rpt 



-A, sf --{stzt~sfzf)\n\zf + s 



B , „B\ 
t ^ 



1 r 
3 



-I i stzt - sf zP + a/^' 



StZf 



2 A. 



where is given in Eqs. (Esl), st, zt, sf , zf , are given in Eqs. (p6|), and 



In U- 



^ ' I 



(60) 



A, = Zt - zf = 



B _ 2Mi [E^ - E) EP 



M2M( 1 



2E I m2 



A, = St - sf = ^ I 1 



M2 



Ml 




(61) 



B. Three body region 

The TBR formula for the electron energy spectrum becomes 

^2 1 f 7 r J / 17 

dr% [E) = ^_MiM2 ^ j A',dz + ^ E 4^*^2'^fe + J H^'Oo + K<^>'^ + E 

where Tj, = 2, r„ = 3, and theej^,* 's are displayed in Eq. (^ ) and in Eq. (|55|). 
The explicit values for i?^, in Eq. (p6|) for fc = 0, 1, 2, 3 are 



si n 



H[6 



"""""^ \ dz 



(62) 



Rn — — 



Az{Tr-2)-As In 



Zt + St 



Zb + \Sb\ 



(63) 



Ri 



(z^ - zl) [Tr - 2) - AsA\s\ - (stzt - SbiZb) In 



Zt + St 



Zb + \Sb\ 



A|s| = St - \sb\ 



(64) 



R2 



zl) {Tr - 2) - 



6 2 



As (zt Sf - Zf, |sfc| 



{stZt - SbiZb) Ai,| + 



stz^ - SbiZ^ ^2^'' 



In 



Zt + St 



Zb + \sb\ 



(65) 
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(zt^-z^) (T;v^-2)-^(zt + z,) 



2 1 

-AsA|s| + -As {stz^ ~ \sb\ zf) + {stz'l - sbizl) Aj^i 



1 



i (stZt - SbiZb) {zt St - Zh \sb\) + (stz^ - Sbizl + i [stZt - SbiZb) ] In 



Zt + St 



Zb + \sb\ 

where Tr is given in Eq. (|57|), and st, zt, Sbi, Zb, A^, and As are given in Eqs. (^6|), and Eqs. (|6l|), respectively. 



(66) 
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VI. CONCLUSIONS 



The singularities contained in the precise analytical formulas at the level of the Dalitz plot are isolated and well 
understood. The logarithmic divergences disappear and the lepton energy spectrum becomes finite when the formulas 
are integrated in an analytical way. The final and general results in Eq. (^9|) and Eq. (^2|) are precise formulas suitable 
to be evaluated numerically, without ambiguity, at any energy in which the charged lepton is emitted. They contain 
the bremsstrahlung in the four body decay region where the three body decay does not take place and also include 
events in which the electron is collinear to the produced hadron (at the edge of the DP). 

Other authors Ref. jl^ have published numerical data for the percent contributions due to the radiative corrections 
to the HSD decay. Only as a test, we have compared booth results. We consider the numerical values obtained with 
the formulas in Eq. (|59| ) and Eq. (|6|) in six Tables. Tables I-III are devoted to the charged HSD process, considering 
the case : S^^n + e + Pe , and Tables IV- VI are entailed to the neutral HSD case, the process we use to illustrate 
this case is A ^ p + e~ + Pe- In the Tables I, HI, IV and VI we exhibit the data given in Ref. for the relative RC 
in %, caused by bremsstrahlung events, which fall inside and outside the TBR Dalitz plot, and the numerical values 
obtained by means of Eq. (^) and Eq. (|6^). The results are in good agreement taking into account the fact that the 
earlier numerical results contain unspecified approximations especially at the Dalitz plot boundary. In addition we 
illustrate our final results with four graphs. For the (S^,n) case, in Fig. |[ the RC at the TBR (which contains the 
virtual and bremsstrahlung corrections) and the RC in FBR (which contains only the bremsstrahlung corrections) 
contributions are displayed separately. In Fig.B, the complete energy spectrum is plotted with and without the total 
radiative corrections. For the (A,p) case. Fig. y^contains the radiative contributions including the Coulomb effect and 
in Fig. ^ we consider the energy spectrum as in the previous case. The Coulomb effect is included in the analysis. It 
is incorporated into the virtual radiative corrections to the energy spectrum for the neutral HSD in the TBR. 

In summary, the analytical results are useful to obtain information derived from experimental data, about the 
underlying interactions in the decay processes, the basic symmetries, and the internal structure of hadrons, through 
the derivation of precise values of the form factors involved in the effective interaction. 

The knowledge of the energy spectrum of the electron is fundamental for the determination of the decay rate in 
these processes. For the complete determination of the decay rate it is necessary to add the contributions of the events 
in the three body and in the four body regions. 

Let us mention that the RC were also computed by other means for photon bremsstrahlung calculations in semilep- 
tonic decays Ref. |lj]. At last, though the evaluation of the radiative corrections for the HSD is a complex and an 
old problem (see list of references in Ref. [Q), the results in Eqs. ( |59| ) and (^2|) are new, and are the culmination of 
a systematical approach to the decays in the whole region of the Dalitz plot. 
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/{, g'l are effective form factors, see Ref. [Q. 
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The FBR model independent functions 6'^^ and 0""^ are given in Eqs. (27) and (p9|), respectively, 
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In this appendix we describe the integration procedure we foUow to obtain the result in Eq. (^6|) , Ref. 
The first step is an integration by parts 
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To take into account the integration limits we consider that 
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+ bozt + = I 



m< E < Ec, 
Ec < E < Em- 



As a consequence of these results iov m< E < Ec and Ec < E < Em we consider two regions. 
We define 



T> 



m < E < Ec, 

Ec < E < Em- 



At the region m < E < Ec taking the limit z ^ Zt 



T§{zt) = 



Zt 



2 (n + 1) z^zt 



lim < In 



(ao + boz + \J z'^ - l)" 



ZtZ-\~ StS 



^n+l 



2 (n + 1) 



[ln|l - ZbZt - SbiSt\ 



,n+l 



2 (n+1) 



In 



2 {boSt + ZtY 



,n+l 



2 (n + 1) 



[ln|l - ZbZt - SbiSt\] ; 



for the limit z ^ Zb 



nizb) 



,n+l 



2 (n + 1) 

Zb 

2 (n+1) 



In 



In 



(ao + boZb + SbiY 



^n+l 



(2s6i) 



-Hi 



^n+1 



2 (n + 1) 



2 (n + 1) 
\n\ztZb - 1 - StSbi\ - 



In \ztZb - 1 - StSbi 



2 (n + 1) 



[ln|l- 
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The following remains: 



-In 



2 (n + 1) 

Similarly for the second region in which < E < 



2 {bost + zt) {ztZb - 1 - StSbi) 



2 (n + 1) 



[In I -2(1 - ZbZt - SbiSt)\\ . 



^n+l 



2 (n + 1) 2^2t 



lim < In 



(oo + 60-2 + - ly 



2(n+l) 



In 



2 (6oSt + Zt) 



ZtZ-1- StS 



2(n+l) 



[ln|l - ^fc^t + SfedStl] 



and 



T^izb) 



2 (n + 1) 



lim < In 



2(n+l) 
ao + boz + Vz'^ - l) 



[ln|l - ZbZt - SbiStW , 



2 (n + 1) 



-In 



1- ZbZ + SbdS 



2 (n + 1) 



[In \ztZb - 1 - StSbd\, 



-2 (6oSfcd + 



2 (n + 1) 



[InjztZfc - 1 - StSfcdl] • 



Then 



T>{zi)-T>{zb) 



Summarizing 



2(n-hl) 



-In 



2 (6oSt + 2;t) {ztZb - 1 - StSftrf) 



2(n+l) 



-In 



(1 - ZbZt - SbiSt) (-2) {boSbd + Zb) 



T^izt)-T<{zb) = 



2 (n + 1) 



A{Em-EY 



J^2 (l _ 2E_ , rn^ 
^*^2 l^-L Ml + 



In general 



where 



T>{zt)-T>{zb) 



n+l ^n+l\ 



2 (n + 1) 



In 



4 - Ef 



+ In 



m2 - eY 



Mi I 



Tn = \n 



Toizt) - Toizb) 



4{Em-EY 



2(n+l) 



■Tr 



2Q{E- Ee) In 



Ml {Em - E) 



Mo 



and <d{E — E^) is the Heaviside function. 
Finally, 



Rn{z) = 



yU+l _ n+1 



1 



2 (n -I- 1) ^ 2 (n -I- 1) 



E 

r=0 



r I ^rA f ^n — r-\-l n— r+1 



- Zf. 



(n-r + l) 



+ {StZ; 



l-SbiZl) j 

J Zh 



-dz. 



FIG. 1. Radiative corrections for the CHSD in the three and four body regions of the Dalitz plot. 



FIG. 2. Electron energy spectrum for the CHSD with all RC contributions and without any. 
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FIG. 3. Radiative corrections for the NHSD in the three and four body regions of the Dalitz plot. 



FIG. 4. Electron energy spectrum for the NHSD with all RC contributions and without any. 



TABLE I. Radiative corrections in the four body region of the Dalitz plot. 



E-(pi) ^ n{p2) + e- (e) + Ue{Pu) + J {k) 


X = E/E^ 
% in Ref. H 
% from Eq. (||) 




0.1 0.2 0.3 
7.8 1.5 0.5 
7.776 1.533 0.458 




0.4 
0.1 
0.135 




0.5 

0.02 

0.019 




TABLE II. Radiative corrections in the three body region 


of the Dalitz 


plot. 








S-(pi) 


n{p2) + e~{£) + Veijpu) 










X = E/Em 

% from Eq. (||) 


0.1 

10.442 


0.2 0.3 0.4 0.5 
5.729 3.329 1.614 0.167 


0.6 

-1.220 


0.7 
2.647 


0.8 

-4.291 


0.9 

-6.597 


TABLE III. Total radiative corrections. 




s-(Pi)- 


n{p2) + e- {£) + u4p^) + J (k) 










X = E/E^ 
% in Ref. jl| 
% from (Eq. Q+Eq. Q) 


0.1 0.2 0.3 0.4 0.5 
18.3 7.1 3.6 1.4 - 0.2 
18.218 7.262 3.787 1.749 0.186 


0.6 
-1.7 
-1.220 


0.7 
-3.2 
-2.647 


0.8 
-5.0 
-4.291 


0.9 
-7.5 
-6.597 




TABLE IV. Radiative corrections in the four body region 


of the Dalitz 


plot. 






A{pi)^p+{p2) + e~(£) + u,(p,)+-y{k) 


X = E/E,^ 
% in Ref. H 
% from Eq. (^) 




0.1 0.2 0.3 
9.5 2.3 0.8 
9.277 2.200 0.765 




0.4 

0.25 

0.242 




0.5 

0.02 

0.024 




TABLE V. Radiative corrections in the three body region 


of the Dalitz 


plot. 








A(pi) - 


P'^ {P2) + e- (e) + u,Xp,^) 










X = E/Em 

% from Eq. (||) 


0.1 

9.058 


0.2 0.3 0.4 0.5 
5.213 3.119 1.531 0.091 


0.6 

-1.396 


0.7 
2.933 


0.8 

-4.692 


0.9 

-7.146 


TABLE VI. Total radiative corrections. 




A(pi) - 


P~^{p2)+e-{e)+u4p^)+j{k) 










X = E/Em 
% in Ref. |l| 
% from (Eq. (|9l)+Eq. (H)) 


0.1 0.2 0.3 0.4 0.5 
18.2 7.1 3.5 1.4 -0.3 
18.335 7.413 3.884 1.773 0.115 


0.6 
-1.8 
-1.396 


0.7 
-3.4 
-2.933 


0.8 
-5.2 
-4.692 


0.9 
-7.7 
-7.146 
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